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Abstract. Quasiclassical boundary conditions for electrochemical potentials at the 
interface between diffusive ferromagnetic and non-magnetic metals are derived for the 
first time. An expression for the boundary resistance accurately accounts for the 
momentum conservation law as well as essential gradients of the chemical potentials. 
Conditions are established at which spin-asymmetry of the boundary resistance has 
positive or negative sign. Dependence of the spin asymmetry and the absolute value 
of the boundary resistance on the exchange splitting of the conduction band opens up 
new possibility to estimate spin polarization of the conduction band of ferromagnetic 
metals. Consistency of the theory is checked on existing experimental data. 
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1. Introduction 

Discovery of giant magnetoresistance (MR) in magnetic multilayers |2j , which consist 
of alternating ferromagnetic metal (F) and normal metal (N) layers, stimulated intensive 
studies of spin-polarized transport through layered magnetic structures jHl IH IH] • In a 
paper [Hj the semiclassical theory of MR in magnetic multilayers has been developed for 
geometry, in which the electric current is perpendicular to the layers (CPP geometry). 
It has been shown that when thickness of the metals in a stack is small compared 
with the spin diffusion length, the magnetoresistance of multilayers can be calculated 
within the two-current series-resistor model [3 IHl l^- In this case MR is expressed 
via the resistance of the F/N interface (boundary resistance), = 2[1 =i= ■jp/N]Rp^j^, 
and the bulk specific resistance of the ferromagnetic layer, = 2[1 =i= Pf]p*f 0- In 
these expressions Jf/n and (3f are parameters of spin asymmetry of the boundary and 
bulk scattering resistances. It follows from works [HI that their linear combination, 
PFiP*FitFi + liR*F/Ni determines sign of MR in layered [Fl/ N / F2l N] x n structures {tF 
is the thickness of the ferromagnetic layers). In works jTTl |^ positive (inverse) MR 
due to negative value of "^f/n in {Fl/Cr / F2/Cr) multilayers has been observed for the 
first time. 

Theoretical calculations of the boundary resistance [121 CH CS] established the 
strong influence of the spin-dependent band structure of ferromagnetic metals on 
magnetoresistance. Using the approach of Ref. |13j and transmission probabilities 
through the F/N interface, calculated earlier from first principles fHl, Stiles and Penn 

obtained from numerical calculations a negative value of "^f/n for the Fe/Cr 
interface, and a positive "^f/n for iron-group ferromagnet - noble metal interfaces. 
However, conditions on the parameters of contacting materials resulting in negative 
or positive values of spin-asymmetry of the boundary resistance were not discussed in 
the above papers. 

In this paper we derive for the first time quasiclassical boundary conditions for 
electrochemical potentials of diffusive ferromagnetic and normal metals, which can be 
used for solution of a wide class of problems in spintronics. Our expression for the 
boundary resistance accurately accounts for the momentum conservation law as well as 
essential gradients of the chemical potentials. We establish conditions on parameters of 
the contacting metals, at which spin asymmetry of the boundary resistance has a positive 
or negative sign. Dependence of the spin asymmetry and the boundary resistance on 
the exchange splitting of the conduction band offers one more way to estimate spin 
polarization of conduction band of ferromagnetic metals. We give an example of such 
an estimation. 

2. Boundary conditions for electrochemical potentials 

We derive boundary conditions for electrochemical potentials of diffusive metals using 
the quasiclassical Green functions technique. The "diffusive ferromagnetic metal" is 
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a metal in which spin sphtting of the conduction band is small compared with the 
momentum relaxation rate for conduction electrons. Let us consider that axis x is 
perpendicular to the F/N boundary, and neglect reversal of the electron spin upon 
transmission through the interface. Then, for each of the metals equations for the 
Green functions ga{'n.,x, p,t) read [T7] : 



dga,a , dgs 

+ Vll 



+ — {9s,a - 9s,a) = 0, (1) 



OX op 
dgs,a , dg 

Here n = Px,a/\Pa\; 9s{a),a = l/2[ga{n^,x,p,t) ± g ^i-n^, x, p,t)] is the single-particle 
quasiclassical Green function symmetric (antisymmetric) with respect to a projection 
of the Fermi momentum „ on the axis x ; is a projection of the Fermi velocity on 
the axis x; a = (|, |) is a spin index, and p = {y,z) is a coordinate in a plane of the 
contact. The bar above gs,a means integration over the solid angle: "gg a — § dVt/2ii gg^^- 
The boundary conditions to Eqs. ((H) and Q follows |T7j: 

9aA^) - 9aA^) - \ 0, min(p^,p^) < p\\ ' 

2Ra.9a,M = Da (^?;.(0) - g^jO)) . (4) 

In the equations Q and (0)) and p'^ are the Fermi momenta in ferromagnetic and 
normal metals, respectively; p\\ is a projection of a momentum on the plane of the 
contact; Da and Ra = l — D^ are the spin-dependent, quantum-mechanical transmission 
and reflection coefficients. Boundary conditions 0, Q obey the specular reflection law: 

p\\ = p^ sinOi = p^ sm6i^ = p^ smO^- (5) 

The angles 6' in (0) are measured from the axis x, a range of variation for the biggest 
one is [0,7r/2]. The quasiclassical equations ((T]) and (j2I), and the boundary conditions 
(121 ) and (jH) are formulated for a single electron trajectory determined by the angles if 
and 6. 

Upon solution of the system of equations (|T)) and Q we shall consider that the 
ferromagnet is located to the left of the boundary x = 0, and the normal metal 
- to the right {x > 0), and that the functions gs,a are homogeneous in the plane 
of the contact. Then, the system of equations (P) and (j21) can be solved in the 
form of integral equations for the functions ga,a and gs,a in the energy representation, 
(?,,„(£) =2 tanh(£/2T) + /,,„(£): 

oo 



fsM = 9a J^) + r- 1 d^^~-fsAO. (6) 



= -<J^) + 7^ / die^-flM- (7) 
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In a dirty metal the solid- angle averaged function fg ^iO obeys the diffusion equation 
with a decay length which is much longer than the mean free path /q. Then, we 
expand /fi^\0 the r.h.s. of Eqs. (jHI) and ((Zj) near point x and take out from 
the integrals independent on ^ terms. Substituting the resulting expansions into the 
boundary condition (jH) we find: 





(8) 



x=0 



To formulate boundary conditions for the functions f^lf^^ (which are, in fact, chemical 

From 



(9) 



potentials - see below) we use a matching procedure proposed in Ref. 
equation ([T)) it follows that for distances of the order of l^^ a from the interface 



dga,c 
dx 



Hence, 



0. 



C = const 



(10) 



''x,a 9a, a 

in each of the metals. Now we calculate, for example, using the expression (jHJ 
for ga,a{x = 0). Then, we calculate far from the interface using an approximate 



expression for ga^ai^)^ 



9a, 



[X 



-I 



dfs.a 



dx 



which follows from equation (0) after expansion of g^^ on Legendre polynomials. 
Equating values of the constant calculated in the two ways, and applying relationship 
between the averaged Green function and the electrochemical potential, « = (2/vr)/XQ, 
we receive the boundary condition for the electrochemical potentials at the interface 
x = 0: 

.dfi^iO) 



where 



dx 



(/^^(0)-/if(0)) 



(12) 



2,a 



X 




COS( 



COS 



cos(M 



The limits of angular integration must satisfy specular reflection conditions at the 
interface, Eq. (jj)). When deriving Eq. ()12|1 we have used conservation of the current 
density at the interface, which follows from Eq. 



(13) 



jf(0) 



a 



N 



dt^^iO) 



e dx e dx 

where and are the bulk, spin-channel conductivities of the metals: 



(J. 



FiN) 



(14) 



(15) 
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Equation (fT^ is actually the second, complementary to Eq. (fT^. boundary condition 
for the semiclassical description of the spin-polarized transport in magnetic multilayers 
in terms of electrochemical potential. 

3. Resistance of the interface 

The derivative from the electrochemical potential in Eq. ()12j) can be expressed in terms 
of the density of current (JU, and we find the spin- dependent resistance of the interface 

/^a(0)-/^f(0)=er,j„, (16) 
^ I- S,,^ 

where A is the area of the contact. It follows from equation (jl7|) that, in the quasiclassical 
approach at specular reflection from the interface, the boundary resistance between 
ferromagnetic and normal metals is determined only by the Fermi momenta of the 
contacting metals and coefficient of transmission through the interface. 

Experimental data are given for the spin asymmetry of boundary resistance, Jf/n, 
and for the renormalized resistance of the interface, ARp^^, determined as follows: 

JF/N = '^^^, AR*p/^ = ^ir^ + r^). (18) 

To calculate dependence of ■Jf/n and AR*p^j^ on the Fermi momentum of the non- 
magnetic metal for various values of the ferromagnet conduction band polarization 
we have used the Fermi-momentum-mismatch model for the transmission coefficient: 
Da = 4f^cff^Q/[(f^„)^ + (f^^a)^]. Results are presented in Figures 1 {p^ > pf > p^), 
2 {p^ > p^ > pf) and 3 {p^ > pf > pf). From our calculations it follows that for 
a non-magnetic metal with a low density of conduction electrons (small value of the 
Fermi momentum p^) the spin asymmetry of the boundary resistance 'Jf/n is always 
negative (Fig. 1). On the contrary, for a non-magnetic metal with a high density of 
conduction electrons the spin asymmetry of the boundary resistance is always positive 
(Fig. 3). In an intermediate situation •Jf/n can be negative as well as positive (see Fig. 
2). To attain the maximum amplitude of negative magnetoresistance in F/N multilayers 
the spin asymmetry of the boundary resistance Jf/n and the asymmetry of the bulk 
resistance Pf should be both positive and close to unity. According to our calculations, 
the Fermi momentum of the non-magnetic metal should be as far as possible close to the 
Fermi momentum of the majority subband of ferromagnetic metal (Fig. 2, p^ /pf 1.0 
and Fig. 3, pf /p^ 1.0). The spin asymmetry of the bulk resistance can be adjusted by 
the type and concentration of impurities in the ferromagnetic metal ^Hl- Clearly, similar 
arguments can be applied to the opposite case of negative values of 'Jf/n and iSf, which 
will result in positive magnetoresistance in multilayers of alternating ferromagnetic and 
non-magnetic metals. However, a negative asymmetry of the bulk resistance is met 
less often ^H]- Competition of opposite in a sign asymmetries of boundary and bulk 
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resistances can result in negative or positive magnetoresistance depending on the choice 
of materials and thickness of the ferromagnetic layers fT^ . 

4. Discussion of experiments 

Experiments on CPP transport in multilayers are very complicated because the 
resistance of a stack of layers of nanometer thickness is very small (order of fQ 
m^). Nevertheless, available experimental data (see reviews jU 120] and references in 
them, and also works f21Ell22]) allow to test the internal consistency of the theory. 
Multilayers of the iron-group ferromagnetic metals with noble metals most likely belong 
to case 3 {p^ > > pf), and parameter of the spin asymmetry 'Jf/n is positive 
jllliai2niiniE21- For example, -fco/Cu ^ 0.77 @]. Then, intersection of the horizontal 
dash line jco/Cu — 0.77 in the top field of Fig. 3 with the curve ■Jf/n{p'^), corresponding 
to 6 = 0.6, gives pf /p'^ — 0.7. Accepting p^ = 1.0A~^ as a trial value for the Fermi 
momentum of the majority subband of cobalt we receive p*"" ~ 1.41^4"^, which fits fairly 
good the free-electron-model value for copper, Pfem — 1-36A~^ jSB]- There are data for 
silver as the non-magnetic spacer: 'Jco/Ag — 0.85 [22] ■ In a similar way, we obtain from 
Fig. 3, p"^^ ~ 1.22y4~^, which fits well the free-electron-model value Pf^em — 1-20A~^ 

Let us look now at consistency of the theory with boundary resistance data. For the 
combination Co/Cu, Ai?^g^^^(exp) ~ 0.51 /fim^ Continuing the vertical dash line 
for copper in Fig. 3 into the bottom field till intersection with the curve AR*p^j^{p^), 
corresponding to 6 = 0.6, we obtain AR^^^(j^{theoT) ~ 0.74 /fim^. A similar procedure 
gives for silver AR'^^^j^g^theoi) ~ 0.69 ffl rn^ (compare with y4i?^^y^^(exp) ~ 0.56 fVL rn^ 
|22| IT2]). It is worthy to notice that the theory reproduces fairly good closeness of the 
boundary resistances of the Co/Cu and Co/Ag interfaces. 

For Co/Cr multilayers the asymmetry of the boundary resistance is negative, 
Ico/Cr — —0.24 [12]. This value is admissible for cases 1 (Fig.l, top field) and 2 
(Fig. 2, top field). Without details, we conclude that the first case does not match 
the expected value of the conduction band polarization parameter for Co 5 ~ 0.6 ±0.1 
as well as results in about two orders in magnitude higher boundary resistance. The 
second case (see dash lines in Fig. 2) results in a Fermi momentum p"-^*" ~ 0.68^4"^, 
and in a boundary resistance Ai?^^^,^^ (theor) ~ 0.61 /i7m^. The Fermi momentum is 
satisfactory in the frame of the free electron model |23^. The boundary resistance is 
close enough to the experimental value y4i?^^yp^(exp) ~ 0.48 /f2m^ [12] ■ We expect 
that better matching of the band structures of cobalt and chromium, both belonging to 
the iron-group metals, would results in a weaker infiuence of the real band structure on 
the boundary resistance. 

A discrepancy with experiment in an absolute value of boundary resistance of about 
20-45% seems not catastrophic because of the following reasons. First, the trial choice 
of p^ = 1.0A~^ was not optimized. Second, we used the free s-electron model as a 
background for the theory. One may expect that for an interface between a metal 
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with predominantly d-electron conduction band (iron group) and an s-electron metal 
(Cu, Ag), reduced overlapping and symmetry mismatch may noticeably increase the 
boundary resistance. Third, Garcia and Stoll have shown (Fig. 2 of Ref. j21]) that the 
interface roughness also increases boundary resistance. This increase is estimated below 
20-60% for different models and magnitudes of the interface roughness, and reasonable 
differences in the Fermi momenta of contacting metals. The interface roughness is 
much less important for the case of GPP transport, considered in this paper, than for 
current-in-plane (GIP) transport [2S1- Finally, spin-reversal at the interface opposes 
the both previously considered processes, decreasing the boundary resistance. The 
reversal of the electron spin by the spin-orbit interaction is always expected upon 
refraction of the electron wave or scattering on roughness at the interface between two 
metals. Quantitative analysis of the competition between the above minor mechanisms 
of boundary resistance is beyond the scope of the paper. However, it is our expectation 
that the key quantity, spin asymmetry of the boundary resistance 'Jf/n, is only weakly 
dependent on band structures matching, interface roughness, spin reversal etc. because 
of considerable cancellations in the dimensionless ratio, Eq. (fTHj) . 

Our trial evaluations show that the experimental data for the spin asymmetry of 
boundary resistance and the absolute value of boundary resistance in Co/Cu, Co/Ag 
and Co/ Cr multilayers can be consistently described with the use of the spin polarization 
parameter for the conduction band of cobalt 5 ~ 0.6. At the level of the experimental 
accuracy and completeness of the theory the estimated value of 5 is identical to 5 ~ 0.57, 
which we have estimated |2E] from experiments of Garcia et al. |2ZI on magnetoresistance 
of cobalt nanocontacts. Experiments on Andreev spectroscopy cite similar values of 5 
for cobalt |2H]. Thus, the spin asymmetry of the boundary resistance in combination 
with the absolute value of the boundary resistance can be used for estimations of the 
spin-polarization parameter 5 of the conduction band of ferromagnetic metals. 

Acknowledgments 

The work was supported by EG via the Gontract No NMP4-GT-2003-505282, and by 
RFBR via projects No 03-02-17432 and No 03-02-17656. L R T acknowledges numerous 
illuminating discussions with Professor Nicolas Garcia. 

References 

[1] Baibich M N , Broto J M , Fert A, Nguyen Van Dau F, Petroff F, Eticnnc P, Creuzet G, Friederich 

A and Chazelas J 1988 Phys. Rev. Lett. 61 2472. 
[2] Binash G, Grunberg P, Saurenbach F and Zinn W 1989 Phys. Rev. B 39 4828. 
[3] Gijs M A M and Bauer G E W 1997 Adv. Phys. 46, 285. 
[4] Bass J and Pratt W P Jr. 1999 J. Magn. Magn. Mater. 200 274. 
[5] Zutic I, Fabian J and Das Sarma S 2004 Rev. Mod. Phys. 76 323. 
[6] Valet T and Fert A 1993 Phys. Rev. B 48 7099. 
[7] Zhang S and Levy P M 1991 J. Appl. Phys. 69 4786. 
[8] Bauer G E W 1992 Phys. Rev. Lett. 69 16676. 



Boundary resistance in magnetic multilayers 



8 



[9] Lee S F, Pratt W P Jr., Yang Z, Holody P, Loloee R, Schroeder P A and Bass J 1993 J. Magn. 
Magn. Mater. 118 LI. 
[10] Rashba E I 2002 Eur. Phys. J. B29 513. 

[11] Vouille C, Pert A, Barthelemy A, Hsu S Y, Loloee R and Schroeder P A 1997 J. Appl. Phys. 81 
4573. 

[12] Vouille C, Barthelemy A, Mpondo F E, Pert A, Schroeder P A, Hsu S Y, Reilly A and Loloee R 

1999 Phys. Rev. B 60 6710. 
[13] Schcp K M, van Hoof J B A N, Kelly P J, Bauer G E W and Inglesfield J E 1997 Phys. Rev. B 

56 10805. 

[14] Stiles M D and Penn D R 2000 Phys. Rev. B 61 3200. 

[15] Xia K, Kelly P J, Bauer G E W, Turek I, Kudrnovsky J and Drchal V 2001 Phys. Rev. B 63, 
64407. 

[16] Stiles M D 1996 J. Appl. Phys. 79, 5805; Phys. Rev. B 54 14679. 
[17] Zaitsev A V 1984 ZhETF 86 1742 [1984 Sov. Phys. JETP 59 1015]. 

[18] Kupriyanov M Yu and Lukichcv V F 1988 ZhETF 94 139 [1988 Sov. Phys. JETP 67 1163]. 
[19] Pert A and Campbell I A 1976 J. Phys. F 6, 849. 
[20] Bass J and Pratt W P Jr. 2002 Physica B 321 1. 

[21] Zambano A, Eid K, Loloee R, Pratt W P Jr. and Bass J 2002 J. Magn. Magn. Mater. 253 51. 
[22] Lee S-F, Yang Q, Holody P, Lolocc R, Hcthcrington J H, Mahmood S, Ikegami B, Vigen K, Henry 

L L, Schoeder P A, Pratt W P Jr. and Bass J 1995 Phys. Rev. 52 15426. 
[23] Ashkroft N W and Mermin N D 1976 Solid State Physics (Saunders, Philadelphia) Table 2.1. 
[24] Garcia N and StoU E 1988 Phys. Rev. B 37, 4415. 
[25] Garcia N, private commimication. 

[26] Tagirov L R, Vodopyanov B P and Efetov K B 2002 Phys. Rev. B 65 214419. 
[27] Tatara G, Zhao Y-W, Munoz M and Garcia N 1999 Phys. Rev. Lett. 83 2030. 
[28] Soulen R J Jr., Osofsky M S, Nadgorny B, Ambrose T, Broussard P, Byers J, Tanaka C T, Nowack 
J, Moodera J S, Laprade G, Barry A and Coey M D 1999 J. Appl. Phys. 85 4589. 

Figure captions 

Fig. 1. Dependence of the spin-asymmetry of boundary resistance '"fp/N (top field) 
and the renormahzed resistance of the interface AR*p^j^ (bottom field) on the Fermi 
momentum of the non-magnetic metal for the case {p^ > > p^). 

Fig. 2. The same as in Figure 1, but for the case {p^ > p^ > p^). 

Fig. 3. The same as in Figure 1, but for the case {p^ > pf > pf )■ 




Case 2: p >p,>p 
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